Future weak lensing constraints in a dark coupled universe 
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Coupled cosmologies can predict values for the cosmological parameters at low redshifts which 
may differ substantially from the parameters values within non-interacting cosmologies. Therefore, 
low redshift probes, as the growth of structure and the dark matter distribution via galaxy and 
weak lensing surveys constitute a unique tool to constrain interacting dark sector models. We focus 
here on weak lensing forecasts from future Euclid and LSST-like surveys combined with the ongoing 
Planck cosmic microwave background experiment. We find that these future data could constrain 
the dimensionless coupling to be smaller than a few xl0~ 2 . The coupling parameter £ is strongly 
degenerate with the cold dark matter energy density Q, c h 2 and the Hubble constant Ho- These 
degeneracies may cause important biases in the cosmological parameter values if in the universe 
there exists an interaction among the dark matter and dark energy sectors. 
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II. COUPLED COSMOLOGIES 
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Cosmological probes indicate that the universe we ob- 
serve today possesses a flat geometry and a mass energy 
density made of ~ 30% baryonic plus cold dark matter 
and 70% dark energy, responsible for the late-time ac- 
celerated expansion [l|-y]- While the ACDM model (a 
fiat universe with a cosmological constant) can describe 
the current observational data, there exist also dynam- 
ical options for dark energy, as the quintessence fluid, 
in which a cosmic scalar field is slowly approaching its 
ground state. This quintessence field, in principle, may 
couple to other fields in nature. Observations strongly 
constrain the couplings to ordinary matter Q. How- 
ever, interactions within the dark sectors, i.e. between 
dark matter and dark energy, are still allowed by ob- 
servations [El-[2lj]. Coupled cosmologies, in order to sat- 
isfy CMB constraints, predict values for the cosmologi- 
cal parameters today which may differ substantially from 
the parameters values within non-interacting cosmolo- 
gies. Therefore interacting cosmologies can hide their 
effects at low redshifts and weak lensing measurements 
can help enormously in constraining dark sector coupled 
models. 

We review here the future CMB constraints on in- 
teracting dark matter-dark energy models presented in 
Ref. 1201 1 , ad ding weak lensing data from the future Eu- 
clid 0, 111] and LSST dj-like surveys. Weak lensing 
probes are shown to be highly complementary to CMB 
measurements and extremely powerful tools to constrain 
interacting dark sector models. 

The structure of the paper is as follows. Section |TT] 
presents the basics of coupled cosmologies. In Sees. IIIII 
and IIVI we describe the lensing extraction methods for 
galaxy surveys and CMB measurements respectively. 
Section [V] contains the description of the future data 
used in our analyses. The results from our Markov Chain 
Monte Carlo (MCMC) analyses are presented in Sec. lVIl 
We draw our conclusions in Sec. IVIII 



At the level of the stress-energy tensor it is always 
possible to introduce an interaction between the fluids of 
the dark sector in the following way [25| : 
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The 4-vector Q v governs the energy-momentum trans- 
fer between the dark components and T^ Am ^ v and Tfcs 
are the energy-momentum tensors for the dark matter 
and dark energy fluids, respectively. The momentum ex- 
change Q v can be parallel to the dark energy four veloc- 
ity uj> dc ^ or to the dark matter four velocity w^ dm ^. The 
first option include all quintessence coupled models and 
are effectively "modified gravity" models, implying the 
presence of a "fifth force" effect (only for the dark mat- 
ter), that is, a violation of the equivalence principle. For 
both options the evolution equations for the dark matter 
and dark energy background energy densities are identi- 
cal and reduce to: 



Pdm + 3%pdm = aQ ! 
Pdc + 3"Hpdo(l + w) = -at 



(2) 
(3) 



where the bars denote background quantities, the dot in- 
dicates derivative with respect to conformal time dr = 
dt/a, % = a/ a and w is the dark-energy equation of 
state. For Q > 0, the energy flows from the dark en- 
ergy system to dark matter one. For Q < 0, the energy 
flow is reversed. In coupled cosmologies the momentum 
exchange can be proportional to the dark matter energy 
density (Q oc pdm) or proportional to the dark energy 
energy density (Q oc pdc)- ■ However, even if models pro- 
portional to the dark matter and dark energy velocities 
provide the same background history, the perturbation 
evolution is dramatically different. Therefore, while ge- 
ometrical probes alone are unable to distinguish among 
the two of them, probes of the perturbation evolution via 
weak lensing measurements will make these two models 
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fundamentally different. Another aspect of coupled mod- 
els is that t hey can show non adiabatic, early time insta- 
bilities [HHi3, E3, EE E3 due to the dark coupling term 
which appears in the dark energy pressure perturbations. 
In the following, we shall restrict our analyses to coupled 
models which satisfy the stability criterion of Ref. [lH 
and therefore are free of early-time, non adiabatic insta- 
bilities. We consider the dark coupled model of Ref. [l4| 
(see also Ref. [18[ for the perturbation analysis details) 



Q 



(dm) 



£Hp dc u 



(dm) 



e) 



(4) 



where £ is a dimensionless coupling (considered constant, 
as well as the dark energy equation of state w, in the 
present analysis). H and pde refer to the total expansion 
rate and dark energy density, background plus perturba- 
tion, i.e H = H/a + 5H and pdo = Pdo + ^Pdo respectively. 
Notice from Eq. that Q„ has been chosen parallel 

to the dark matter four velocity u„ , in order to avoid 
momentum transfer in the rest frame of the dark matter 
component [ll|. For this choice of energy exchange Q\?\ 
positive (negative) values of the coupling £ will lead to 
lower (higher) dark matter energy densities in the past 
than in the uncoupled £ = case. We only consider here 
negative couplings and w > — 1, avoiding the instability 
problems previously mentioned, see Ref. [l4j for details. 
For the numerical analyses presented here, we have mod- 
ified the publicly available CAMB code [2g|, taking into 
account the presence of the dark coupling in both the 
background and the linear perturbation equations. 



III. GALAXY WEAK LENSING 

Weak gravitational lensing of the images of distant 
galaxies offers a useful geometrical way to map the 
matter distribution in the Universe. Following Ref. [29[ 
one can describe the distortion of the images of distant 
galaxies through the tensor: 
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where K is the convergence term and 7 = 71 + 172 is 
the complex shear field. As shown in Ref. [3(| the shear 
and the convergence terms can be written as a function 
of the projected Newtonian potentials ipjj". 
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where the commas indicate the derivatives with respect 
to the directions transverse to the line of sight and 
the projected potentials are ip^j = —(1/2) J g(z)("9 + 
<S>^ij)dz with the lensing kernel : 

,n{z!)D A (z,z!) 



9{z) 



dz' 



Here n(z) is the galaxy redshift distribution. In our anal- 
ysis we assume flatness of the Universe. However in gen- 
eral the angular diameter distance Da between the lens 
and the source depends on the spatial curvature K: 

D A = — L sin(\/Fr), K > 
v K 



K = 



D A = r, 



D A = 



1 



i—K 

and the comoving distance is 

r(z, z ) 



sinh( v /Z ^r), K < 



dz' 

m7] 



with E{z) = H(z)/H a . 

Images distortions induced by the matter distribution 
are generally small. To extract cosmological information 
it is hence necessary to statistically analyze a large num- 
ber of images. The two point correlation function of the 
convergence is at present the best measured statistic of 
the weak lensing but, of course, also higher order statis- 
tics contains cosmological information. It is convenient 
to work in the multipole space and define the conver- 
gence power spectrum as the harmonic transform of the 
two-point correlation function. This is usually the most 
analyzed and studied statistical quantity related to the 
weak lensing and we will focus on the convergence power 
spectra in order to properly compare our results to sim- 
ilar analysis in literature. However it should be stressed 
that, as shown in 31|, the convergence power spectrum 



is only indirectly and partially obtainable from the two 
point correlation function. 

Future surveys will measure redshifts of billions of 
galaxies allowing the possibility of a tomographic recon- 
struction of the matter distribution. We can define hence 
the convergence power spectra in each redshift bin and 
the cross-power spectra: 



f°° dz 

PjkW - H$ ^ -^ Wt (z)W J (z)P NL [ 



r(z) 
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where Pnl is the non-linear matter power spectrum at 
redshift z, obtained correcting the linear one Pl- W(z) 
is a weighting function: 

Wi{z) = |0 m (l + z) f Z ' +1 dz' ^^ (6) 
2 Jz z r(0,z') 

with subscripts i and j indicating the redshift bin. Equa- 
tion ([5]) shows the cosmological information contained in 
weak lensing measurements: the function W(z) encodes 
the information on how the three-dimensional matter dis- 
tribution is projected on the sky, while the matter power 
spectrum quantifies the overall matter distribution. 
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The observed convergence power spectra is affected 
mainly by systematic uncertainties arising from the in- 
trinsic ellipticity of galaxies "f^ ms - These uncertainties 
can be reduced averaging over a large number of sources. 
The observed convergence power spectra will be hence: 
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(7) 



where hj is the number of sources per steradian in the 
j — th bin. 



IV. CMB LENSING EXTRACTION 

The analysis presented here includes, in addition to 
the primary CMB anisotropy angular power spectrum, 
the information from CMB lensing. Gravitational CMB 
lensing, as already shown (see e.g. [HI, [33|) can improve 
significantly the CMB constraints on several cosmolog- 
ical parameters, since it is strongly connected with the 
growth of perturbations and gravitational potentials at 
redshifts z < 1 and therefore, it can break important de- 
generacies. The lensing deflection field d can be related 
to the lensing potential as d — V<p [34[. In harmonic 
space, the deflection and lensing potential multipoles fol- 
low 



and therefore, the power spectra C\ 
df* = ((f>T4>?*) are related through 



dd — 



(8) 
) and 



V. FUTURE DATA ANALYSIS 
A. Galaxy weak lensing data 

Future weak lensing surveys will measure photometric 
redshifts of billions of galaxies allowing the possibility of 
3D weak lensing analysis ie.a.l^Tt ■'>!) : or a tomographic 
reconstruction of growth of structure as a function of time 
through a binning of the rcdshift distribution of galax- 
ies, with a considerable gain of cosmological information 
(e.g. on neutrinos [40J; dark energy [39(; the growth of 
structure jU HH and the dark matter distribution as a 
function of redshift (43J). 

Here we use the typical specifications for future weak 
lensing surveys like those of the Euclid and LSST ex- 
periments. Euclid will observe about 35 galaxies per 
square arcminute in the redshift range 0.5 < z < 2 
with an uncertainty of about <r z = 0.05(1 + z) (see 
LSST is expected to have similar characteristics, with 
slightly higher number density of sources and larger red- 
shift range, but also with a higher intrinsic shear. We 
build mock datasets of convergence power spectra for 
these two surveys. Tables U and QT] show the number of 
galaxies per arcminute -2 (n ga i), redshift range, f s ky and 
intrinsic ellipticity for these surveys. 
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0.5 < z < 2 0.5 0.22 



cr = i(i + i)cf 



(9) 



Lensing introduces a correlation between different CMB 
multipoles (that otherwise would be fully uncorrclatcd) 
through the relation 



TABLE I. Specifications for the Euclid like survey considered 
in this paper. The table shows the number of galaxies per 
square arcminute (n ga l), redshift range, f s k y and intrinsic el- 
lipticity (7rms) P er component. 
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where a and b are the T, E, B modes and H is a linear 
combination of the unlensed power spectra Cf h (see 35 1 
for details). 

In order to obtain the deflection power spectrum from 
the observed Cf b , we have to invert Eq. (fTUl) . defining a 
quadratic estimator for the deflection field given by 



d(a, b) 



, ab 
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where n°£ is a normalization factor needed to construct 
an unbiased estimator (d(a, b) must satisfy Eq. ©). The 
variance of this estimator reads 



(d{a,b)¥'d{a',V)$') 
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(12) 



and depends on the choice of the weighting factor W 
and leads to a noise N£ a hb on the deflection power 
spectrum Cff obtained through this method. In the 
next section we describe the method followed here to 
extract the lensing noise. 



rigai (arcmin 



redshift f sky ^ r 



40 



< z < 3 0.5 0.28 



TABLE II. Specifications for the LSST like survey considered 
in this paper. The table shows the number of galaxies per 
square arcminute (n ga l), redshift range, f s k y and intrinsic el- 
lipticity dispersion ("frms) per component. 



The expected la uncertainty on the convergence power 
spectra P(£) is given by [44| : 



<7( 



(21 + l)f sky A e 



P(l) 



7 2 
Irms 

n ga l 



(13) 



where Ai is the ^-bin width used to generate data. 
Here we choose Ai = 1 for the range 2 < I < 100 and 
= 40 for 100 < I < 1500. For the convergence power 



spectra we use £„ 



1500 in order to exclude the 



scales where the non-linear growth of structure is more 
relevant and the shape of the non-linear matter power 
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spectra is, as a consequence, more uncertain (see [45j). 
We describe the galaxy distribution of Euclid survey 
as in 46], n(z) oc z 2 exp(— (z/zo) 1 ' 5 ) where zq is set 
by the median redshift of the sources, zq — z m /\A\. 
Here we calculate the power spectra assuming a median 
redshift z m = 1. Although this assumption is reasonable 
for the Euclid survey, it is known that the parameters 
that control the shape of the distribution function may 
have strong degeneracies with some cosmological param- 
eters as the matter density, cr 8 an d the spectral index . 



We compute the minimum variance lens- 
ing noise for the Planck experiment by 
means of a routine publicly available at 
http : //lesgourg . web . cern . ch/lesgourg/ codes . html. 
The datasets (which include the lensing deflection power 
spectrum) are analyzed with a full-sky exact likelihood 
routine available at the same URL. 



C. Analysis method 



B. CMB data 

We create a full mock CMB dataset (temperature, El- 
polarization mode and lensing deflection field) with noise 
properties consistent with the Planck [48| experiment 
(see Tab. IIIII for specifications). 

Experiment Channel FWHM AT /T 
Planck 70 14' 4.7 

100 10' 2.5 
143 7.1' 2.2 

fsky = 0-85 



TABLE III. Planck experimental specifications. Channel 
frequency is given in GHz, FWHM (Full- Width at Half- 
Maximum) in arc-minutes, and the temperature sensitivity 
per pixel in jj,K/K. The polarization sensitivity is AE/E = 
AB/B = V2AT/T. 



We consider for each channel a detector noise of 
ur 1 = (6a) 2 , where 6 is the FWHM (Full- Width at Half- 
Maximum) of the beam assuming a Gaussian profile and 
a is the temperature sensitivity AT (see Tab. IIIII for the 
polarization sensitivity). We therefore add to each Ce 
fiducial spectra a noise spectrum given by: 



N t = w- 1 exp(t(e + l)/£ 2 b ) . 



(14) 



where If, is given by i b = V8hi2/#. 

We make use of the method presented in [35[ to con- 
struct the weighting factor W of Eq. (jTTJ) . In that pa- 
per, the authors choose VP to be a function of the power 
spectra Cf b , which include both CMB lensing and pri- 
mary anisotropy contributions. This choice leads to five 
quadratic estimators, with ab = TT, TE, EE, EB, TB; 
the BB case is excluded because the method of Ref. [35[ 
is only valid when the lensing contribution is negligible 
compared to the primary anisotropy, assumption that 
fails for the B modes in the case of Planck. 
The five quadratic estimators can be combined into a 
minimum variance estimator which provides the noise on 
the deflection field power spectrum Cf d : 



dd 
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We perform two different analyses. First, we compute 
the expected constraints on the coupling parameter £ 
from Planck and Euclid data, comparing the results with 
the limits arising from Planck and LSST data. Secondly, 
we investigate the effects of a wrong assumption about 
the interaction between dark matter and dark energy on 
the values of the cosmological parameters: we generate 
a dataset with a non-zero £ fiducial value but analyze 
the data assuming that there is no coupling between the 
dark components (£ = 0). We perform a MCMC analy- 
sis based on the publicly available package cosmomc [49[ 
with a convergence diagnostic using the Gelman and Ru- 
bin statistics. 

We sample the following set of cosmological parame- 
ters, adopting flat priors on them: the baryon and cold 
dark matter densities ftbh 2 and ft c h 2 , the ratio of the 
sound horizon to the angular diameter distance at de- 
coupling S , the scalar spectral index n s , the overall nor- 
malization of the spectrum A s at k — 0.002 Mpc -1 , the 
optical depth to reionization r, and, finally, the coupling 
parameter £. 

The fiducial model for the standard cosmological pa- 
rameters is the best-fit from the WMAP seven year data 
analysis of Ref. with n b h 2 = 0.02258, n c h 2 '= 0.1109, 
n s = 0.963, t = 0.088, A s = 2.43 x 10" 9 and 9 = 1.0388. 
For the coupling parameter we first assume a fiducial 
value £ = to test the constraints achievable on the 
coupling model. Finally, we analyse a dataset with a 
fiducial value £ = —0.1 assuming (wrongly) a ACDM 
scenario with £ = 0, with Planck and Euclid forecasted 
data. This exercise will allow us to investigate the bias 
introduced on cosmological parameter inference from a 
wrong assumption about the coupling model. 



VI. RESULTS 

In Tabic |IV] we show the MCMC constraints at 68% 
c.l. for the coupled universe from Planck data alone and 
from Planck data combined with Euclid data. For this 
last case we also fit the data fixing £ to 0, thus perform- 
ing a standard analysis in a universe where dark matter 
and dark energy are not interacting, in order to show 
the importance of the degeneracies introduced by the 
presence of a coupling £ on the other cosmological pa- 
rameters errors. There is a very high level of correlation 
among the dimensionless coupling £ and the parameters 
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Hq and fi c /i 2 in the Planck analysis (see also Figs. [3J and 
13]). When Planck and Euclid data are combined, the de- 
generacy between £ and Hq is broken, leading to a much 
better constrain on the coupling parameter £ than when 
using CMB data alone[20j], as one can notice from Table 
IIVI and Fig. [2] However, the degeneracy between £ and 
VL c h? is not broken by the combination of Planck and 
Euclid, thus it will be possible to further improve the 
constraints on £ with independent measurements of f2 c . 
We also note that the constraints on the standard cos- 
mological parameters are in good agreement with those 
reported in [23j . 





Planck 


Planck+Euclid 


Model 


Varying £ 


C = o 


Varying £ 


£ = o 


Parameter 












0.00013 


0.00013 


0.00010 


0.00010 


A(n c h 2 ) 


0.0299 


0.0010 


0.0024 


0.00055 


A(0.) 


0.0023 


0.00026 


0.00027 


0.00023 


A(r) 


0.0042 


0.0042 


0.0026 


0.0026 


A(n s ) 


0.0030 


0.0031 


0.0029 


0.0027 


A(log[10 10 A]) 


0.013 


0.013 


0.0097 


0.0093 


A(Ho) 


2.28 


0.43 


0.29 


0.27 


A(fi A ) 


0.0614 


0.0050 


0.0062 


0.0026 




> -0.56 




> -0.05 





TABLE IV. 68% c.l. errors on cosmological parameters. Up- 
per limits on £ are 95% c.l. constraints. 



As stated above, future surveys like Euclid will be able 
to tomographically reconstruct the matter distribution. 
Exploiting this possibility would improve the constraints, 
but, as already pointed out in the non tomographic 
analysis can be thought as a conservative estimation of 
the constraints as we are not including systematic effects. 

Table |V] contains both the results from the com- 
bination of Planck and Euclid data and those from 
the combination of Planck and LSST. Notice that the 
results are quite similar. However, the slightly better 
constraints on £ from the Planck plus Euclid combi- 
nation leads to a better measurement of the cold dark 
matter content of the universe than the one performed 
by Planck plus LSST data (see also Fig. [TJ. 
Moreover, the difference between Planck+Euclid and 
Planck+LSST results would be bigger if systematic 
effects are included, as LSST, being a ground based 
survey, will be more affected by these. 

In addition, we have also (wrongly) fitted a mock 
dataset with £ = —0.1 to a non interacting cosmology in 
which the dimensionless coupling vanishes (£ = 0). From 
this exercise we find a consistent bias in the recovered 
best fit value of the cosmological parameters due to the 
strong degeneracies among £ and both the Hubble con- 
stant Hq and the matter energy density f2 m parameters, 
see Tab lVIl Note, from the results depicted in Figs. [31 0] 
and [5] and also from the results in Tab. IVII that the shift 
in the best fit values is, as expected, along the direction 





Planck+Euclid 


Planck+LSST 


Parameter 








0.00010 


0.00010 


A(n c h 2 ) 


0.0024 


0.0026 


A(0.) 


0.00027 


0.00028 


A(r) 


0.0026 


0.0027 


A(n s ) 


0.0029 


0.0029 


A(log[10 10 A s ]) 


0.0097 


0.010 


A(flb) 


0.29 


0.29 


A(fi A ) 


0.0062 


0.0065 




> -0.04 


> -0.06 



TABLE V. 68% c.l. errors on cosmological parameters. Upper 
limits on £ are 95% c.l. constraints. 




0.23 



-0.1 -0.08 -0.06 -0.04 -0.02 

? 

FIG. 1. 2-D constraints on £ and Q m using Planck and LSST 
data (blue contours) and Planck and Euclid data (red con- 
tours). 



of the degeneracy of £ with these parameters. These 
results show that even for a small value of £, the best 
fit values recovered by wrongly assuming that t here 
is no dark coupling are more than 68% c.l. (for some 
parameters at more than 95% c.l.) away from the correct 
fiducial values, and may induce an underestimation of 
both Hq and as and an overestimation of Q c h 2 . In the 
last column in Tab. I VII we show the difference between 
the wrong value estimated fixing £ = and the fiducial 
value, relative to the la error: as expected the largest 
shifts are in the parameters that are directly involved 
in determining the energy momentum transfer between 
components, namely matter and dark energy densities 
and Hubble parameter (see section |H|) . We note that 
also other parameters, as as and n s , have significant 
shifts. 

We conclude, hence, that future analyses of high pre- 
cision data from Euclid and Planck need to consider 
possible deviations from the minimal ACDM scenario 
in order to avoid biases in the measurements of the 
cosmological parameters. 
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U c h 2 


0.1245 ±0.00061 


0.1083 ±0.0024 


0.1109 


T 


0.086 ± 0.0029 


0.090 ± 0.0033 


0.088 


n s 


0.955 ± 0.0014 


0.961 ± 0.0028 


0.963 


H 


69.6 ±0.23 


70.9 ±0.30 


71.0 


Q A 


0.697 ±0.0031 


0.739 ± 0.0066 


0.735 


°"8 


0.752 ± 0.00044 


0.841 ±0.019 


0.82 



TABLE VI. Best fit values and 68% c.l. errors on cosmological 
parameters for the case in which a fiducial model with f = 
—0.1 is fitted to a ACDM model where £ = is assumed. 






FIG. 2. 2-D constraints on £ and fl m using Planck (blue 
contours) and Planck plus Euclid (red contours). 



VII. CONCLUSIONS 
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FIG. 3. Top panel: 68% and 95% c.l. contours in the (Q, m , 
£) plane from Planck data only. Bottom panel: same as in 
the top panel, but for the combination of Planck plus Euclid 
data (note the different scale for the x-axis). 



The current accelerated expansion of the universe is 
driven by the so-called dark energy. This negative pres- 
sure component could be interpreted as the vacuum en- 
ergy density, or as a cosmic, dynamical scalar field. If 
a cosmic quintessence field is present in nature, it may 
couple to the other fields in nature. While the couplings 
of the quintessence field to ordinary matter are severely 
constrained, an energy exchange among the dark matter 
and dark energy sectors is allowed by current observa- 
tions. Interacting cosmologies, in order to fit Cosmic Mi- 
crowave Background observations, predict non standard 
values for the low redshift universe observables. There- 
fore, measurements of the growth of structure and of the 
dark matter distribution via galaxy and weak lensing sur- 
veys offer a unique window to study the nature of the 
dark sectors. 

The major goals of the on-going Planck experiment 
and the next generation of galaxy surveys are to de- 
termine the nature of the dark energy component and 
to measure the remaining cosmological parameters with 
unprecedented precision. In this study we have exploited 



the capabilities of the former experiments to improve cur- 
rent constraints on the dimensionless coupling £, which 
drives the energy flow between the dark energy and dark 
matter sectors [14 ] . 

We have generated mock data for both the Planck ex- 
periment and the Euclid and LSST-like weak lensing sur- 
veys. CMB gravitational lensing extraction has also been 
included in the analysis. For the weak lensing surveys, 
the observable we have exploited here is the convergence 
power spectrum. The mock data have then been an- 
alyzed using MCMC techniques to compute the errors 
on the several cosmological parameters considered here. 
The constraints on the dimensionless coupling parameter 
of the interacting model fully explored here, £, arc greatly 
improved with respect to previous analyses in which only 
CMB probes were considered. We find £ > —0.04 at the 
95% c.l. for the combination of Planck survey and Eu- 
clid weak lensing data. Data from a LSST-like survey 
combined with data from the Planck experiment could 
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FIG. 4. Top panel: 68% and 95% c.l. contours in the (H , £) 
plane from Planck data only. Bottom panel: same as in the 
top panel, but for the combination of Planck plus Euclid data 
(note the different scale for the x-axis). 



as possible scenarios when fitting their data. Finally we 
conclude noting that, as recently outlined by [5l|, cou- 
pled quintessence is expected to have a distinctive impact 
on the formation of structures at galaxy clusters scales. 
In our work we limited our analysis to linear scales, but 
for the analysis of real future data, it will be necessary a 
precise modeling of the effects of coupled quintessence on 
non-linear scales. At the same time, this small scale effect 
will offer a way to constrain coupled cosmologies through 
measurements of the galaxy clusters mass distribution. 
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FIG. 5. 68% and 95% confidence levels from Planck plus Eu- 
clid data when a fiducial cosmology with £ = —0.1 is fitted to 
a non interacting cosmology with £ fixed to (blue countours) 
or to an interacting cosmology in which £ is allowed to vary. 
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